Three-Dimensional Spin Fluctuations in Na .75CoO2 



O 

o 

(N 

> 

O 

ON 



L. M. Helme/'B A. T. Boothroyd, 1 R. Coldea, 1 D. Prabhakaran, 1 D. A. Tennant, 2 A. Hiess, 3 and J. Kulda 3 

1 Department of Physics, University of Oxford, Oxford, 0X1 3PU, United Kingdom 
2 School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife KY16 9SS, United Kingdom 



' Institut Laue-Langevin, BP 156, 38042 Grenoble Cedex 9, France 
(Dated: February 2, 2008) 

We report polarized- and unpolarized-neutron scattering measurements of magnetic excitations 
in single-crystal Nao.7sCo02. The data confirm ferromagnetic correlations within the cobalt-layers 
and reveal antiferromagnetic correlations perpendicular to the layers, consistent with an A-type 
antiferromagnetic ordering. The magnetic modes propagating perpendicular to the layers are sharp, 
and reach a maximum energy of ~12meV. From a minimal spin wave model, containing only 
nearest-neighbour Heisenberg exchange interactions, we estimate the inter- and intra-layer exchange 
constants to be 12.2 ± 0.5 meV and —6 ± 2meV, respectively. We conclude that the magnetic 
fluctuations in Nao.7sCo02 are highly three-dimensional. 

PACS numbers: 75.40.Gb, 74.25.Ha, 74.20.Mn, 78.70.Nx 
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Its potential as a battery electrode material |lj and in 
thermoelectric devices [2J, and most recently the discov- 
ery of superconductivity after hydration have made 
the layered cobaltite Na x Co02 the subject of intense re- 
search in the last few years. The structure comprises 
two-dimensional layers of C0O2 separated by layers of 
sodium ions, with the Co atoms forming a triangular lat- 
tice 0- The superconducting compound Na^CoCvyH^O 
(x « 0.3, y ~ 1.3, T c rs 4.5 K) is of particular interest as 
the first cobalt-oxide based superconductor ||. The lay- 
ered structure and existence of superconductivity over a 
narrow range of doping || near to a Mott insulator invite 
comparisons with the copper-oxide superconductors, and 
evidence is mounting for an unconventional mechanism 
of superconductivity Q . 

In common with the superconducting cuprates, the 
properties of the normal metallic state of Na^OoC^ are 
unusual and show effects due to strong electronic cor- 
relations. One of the outstanding puzzles is the nature 
of the magnetic interactions, which may play a central 
role in the formation of the superconducting state. In 
the range x ~ 0.7-0.95 the susceptibility of Na x Co02 
shows a sharp magnetic transition at T m w 22 K 
Muon-spin rotation (/iSR) measurements confirmed the 
existence of static magnetic order below T m , and placed 
an upper limit of 0.2 /ie on the size of the ordered mo- 
ment |9j . For T > T m the susceptibility can be fitted to a 
Curie- Weiss law plus a constant term, which indicates a 
degree of local character to the magnetism 7] . Such fits 
give negative Weiss temperatures, implying dominant an- 
tiferromagnetic correlations. However, a recent neutron 
scattering study of Nan.7sCo02 found strong ferromae 
netic spin correlations within the cobalt-oxide laye i 
consistent with several theoretical predictions [llj . Until 
now there have been no measurements of magnetic fluc- 
tuations perpendicular to the layers, and the nature of 
the magnetic order below T m has not been established. 

Here we investigate further the magnetic correlations 
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in Nao.75Co02, using both polarized- and unpolarized- 
neutron scattering. The new data reveal strong antiferro- 
magnetic correlations perpendicular to the cobalt-oxide 
layers, consistent with an A-type antiferromagnetic or- 
dering. Surprisingly, the inter-layer exchange coupling 
is found to be similar in strength to the intra-layer cou- 
pling, so that despite the layered structure the magnetic 
interactions are highly three-dimensional. 

The measurements were made on a single crystal of 
Nan 7RC0O2 grown in Oxford by the floating-zone method 
[l2J . For the neutron studies we cleaved a crystal of mass 
~1.5 g from a zone- melted rod. Smaller crystals from the 
same rod were examined by x-ray diffraction, magnetom- 
etry and electron microscopy. The analysis revealed the 
presence of small inclusions of cobalt oxides (CoO and 
C03O4) consistent with previous reports for melt-grown 
crystals [13J. These impurity phases, which amounted to 
a few per cent of the total, were found by neutron diffrac- 
tion to grow epitaxially on the host lattice. Once the ori- 
entation of the impurity crystallites had been determined 
it was straightforward to distinguish the impurity signal 
from that of the host. The anisotropic magnetic suscepti- 
bility of the crystals exhibited an anomaly at T m « 22 K 
for fields applied parallel to the c axis, consistent with 
the magnetic transition observed previously • 

Unpolarized- and polarized-neutron scattering mea- 
surements were performed on the thermal triple-axis 
spectrometers IN8 and IN20, respectively, at the Insti- 
tut Laue-Langevin. On IN8 we employed a Si (111) 
monochromator and a pyrolytic graphite (002) analyser, 
and worked with a fixed final energy Ef — 14.7 meV. To 
increase the count rate both monochromator and anal- 
yser were curved horizontally and vertically for optimum 
focussing. For the polarized-neutron measurements on 
IN20 we used curved Heusler (111) as both monochro- 
mator and analyser, and E{ = 34.8 meV. On both in- 
struments a graphite filter was placed in the scattered 
beam to suppress higher order harmonics. The crystal 
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FIG. 2: (Color online) Left: The magnon dispersion of 
Nao.75Co02 parallel to (00/) measured on IN8 at a temper- 
ature of 1.5 K. Filled circles are from constant-energy scans, 
e.g. Fig.Q^a), and open squares are from constant-/ scans, e.g. 
Fig.^c). The solid curve is calculated from the spin wave dis- 
persion Eqs. [5] and [3] with exchange constant J c = 12.2 meV. 
Dotted lines show the zone boundaries. Right: The magnetic 
structure on which the spin wave model is based, showing the 
two exchange constants J a b and J c , and the spin direction. 



FIG. 1: (Color online) Magnetic scattering from Nao.7sCo02. 
(a) and (b) show scans parallel to (00/) made at a fixed energy 
of 7meV. (c) displays energy scans made at the zone bound- 
ary / = 2.5 at temperatures of 1.5 K and 70 K. The data in (a) 
and (c) were obtained with unpolarized neutrons on IN8 using 
a fixed final energy of Ef — 14.7 meV. The data in (b) are the 
spin-flip (SF) and non-spin-flip (NSF) scattering collected on 
IN20 with the neutron polarization P parallel to the scatter- 
ing vector Q, and with fixed Ef = 34.8 meV. The horizontal 
bars in (a) and (b) indicate the experimental resolution. 



was mounted in a standard helium cryostat, and most 
of the measurements were made within the (100)-(001) 
scattering plane. Here, the notation (hkl) represents the 
wavevector Q = ha* + kb* + lc*, where a*, b* and c* 
are the basis vectors of the hexagonal r ecip rocal lattice 
(the angle between a* and b* is 60 deg) [l4( . 

Our previous studies lij showed that the magnetic 
excitation spectrum of Nao.75Co02 is centred on the T- 
point of the two-dimensional Brillouin zone, i.e. (00), 
corresponding to ferromagnetic correlations within the 
cobalt-oxide layers. In the present work we concentrated 
on the out-of-plane wavevector component of the mag- 
netic fluctuations. 

FigureOJa) shows an example scan parallel to the (00/) 
direction performed on IN8 at a fixed energy transfer of 
7meV. Two peaks can be seen symmetrically either side 
of I = 3. Figure Hfb) displays the same scan but this 
time performed on IN20 with the neutron polarization 
maintained parallel to the scattering vector during the 
scan. In this configuration the spin-flip (SF) scattering is 
purely magnetic, and the non-spin-flip scattering is non- 



magnetic. The two peaks are clearly present in the SF 
channel and absent from the NSF channel. The peaks 
are essentially resolution-limited, as indicated, but are 
less well resolved in Fig. QJb) than in Fig. ^a) because 
of the larger neutron energy used on IN20. We conclude 
that the peaks arise from magnetic excitations. 

The scan shown in Fig. ^a) was repeated for different 
fixed energies between 3 meV and 10 meV. Each scan 
contained two peaks symmetric about (003), with the 
peak separation increasing with increasing energy. In 
addition, energy scans were made at several fixed points 
along the line (00/)). Figure QJc) shows one such scan, 
made at / = 2.5, the zone boundary in the out-of-plane 
direction. The scan was performed at 1.5 K and then 
repeated at 70 K. The prominent peak at ~12meV in 
the low temperature scan has disappeared by 70 K. This 
again confirms the magnetic origin of the scattering since 
magnetic correlations are destroyed with increasing tem- 
perature. 

By fitting Gaussian functions to the peaks in both 
types of scan we constructed the magnon dispersion rela- 
tion. This is displayed in Fig.El There is clearly a mode 
dispersing from (003) with a maximum energy of approxi- 
mately 12 meV. The crystal structure of Na^OoO^ is such 
that no structural Bragg peaks are allowed for positions 
(00/) with odd /. As expected, therefore, no structural 
Bragg peak was observed at (003), but scans made at 
different temperatures revealed no magnetic Bragg peak 
at this point either. 

The simplest spin arrangement consistent with the ob- 
servations is the A- type antiferromagnet shown in Fig. [21 
in which the spins are ordered ferromagnetically within 
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the layers and the layers are coupled antiferromagneti- 
cally along the c axis. Each cobalt ion is taken to have 
the same spin. For no magnetic Bragg peak to appear 
at (003) the spins must be parallel or antiparallel to the 
c-axis. This is because neutrons do not couple to spin 
components parallel to the scattering vector. 

To analyze the three-dimensional dispersion in more 
detail we compare the experimental results with a spin 
wave model containing the minimum number of exchange 
parameters. The Hamiltonian is 



T~L — J a b ^ ' Si • Si' + J c S ' Si ■ Sj 



(1) 



(i,i>) 



where J a (, and J c are intra- and inter-layer exchange 
parameters, respectively, as indicated in Fig. [21 Only 
nearest-neighbour interactions are included in the sum- 
mations, and (i, i'} and denote spin pairs within the 
same layer and on adjacent layers, respectively. 

Standard methods were used to derive the spin-wave 
dispersion and scattering cross-section from the Hamilto- 
nian. The resulting expression for the spin-wave energy 
dispersion is 

Tiwk = 25 \J (Jk - Jk m )[(^k-k m + Jk+k m )/2 - Jk m ], 

(2) 

where 5 is the spin (here assumed to be S = 1/2), k is 
the magnon wavevector, and k m = (001) is the propaga- 
tion vector of the magnetic structure. Jk is the Fourier 
transform of the exchange couplings, given by 

Jk = J c cos(7ri)+J a 6[cos(27r/i)+cos(27rfc)+cos(27r(/i+fe))], 

(3) 

with Jk m , Jk-k m and Jk+k m defined in a similar manner. 

The dispersion relation along the (00/) direction does 
not depend on J ao , so by comparing the spin- wave dis- 
persion to the data in Fig. [2 we can immediately obtain 
a value for J c . The best fit is shown by the solid curve 
on Fig. |21 which is calculated with J c = 12.2 mcV. At 
low energies the data points lie systematically above the 
fitted curve, suggesting the presence of a small gap of 
1-2 meV. Apart from this, the model provides a good 
description of the data. 

The measurements described so far probe only the 
inter-plane correlations. To gain quantitative informa- 
tion on correlations within the planes we apply the model 
to the results obtained previously with a crystal of the 
same composition using the MAPS time-of-flight spec- 
trometer at ISIS [lOj. Figure [3f a) reproduces part of a 
slice through the MAPS data from Ref. [lOj, in which the 
data were averaged over energy transfers E = 8 — 12 meV. 
The configuration used on MAPS is such that the out- 
of-plane wavevector component / varies with E. For this 
slice lw 1, 

To continue the analysis, a MAPS-style data set was 
simulated from the model to allow direct comparison with 
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FIG. 3: (Color online) (a) Neutron inelastic scattering from 
Nao.75Co02 recorded on the MAPS spectrometer at ISIS with 
an incident energy of 60meV 10,]. The map contains data 
averaged over energy transfers of 8-12 meV, and is projected 
onto the (h, k) reciprocal lattice plane of the crystal, (b) 
Simulated intensity using the model described in the text with 
Jab = — 6meV and J c = 12.2 meV. The axis labels correspond 
to the hexagonal reciprocal axes drawn in the figure. 
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FIG. 4: (Color online) Constant energy cuts taken along the 
line marked X in Fig. Open circles show neutron data 
points, while the solid lines are from the simulation. The 
lOmeV and 14meV data have been displaced vertically by 
3 and 6 units, respectively. The horizontal bar indicates the 
resolution. 



the MAPS data. The MAPS data set is an intensity array 
in (Q, E) space, so for each data point in this space the 
simulated intensity was calculated, including the mag- 
netic form-factor and orientation factor . For the cal- 
culation, J c was fixed to the value 12.2 meV determined 
from the inter-layer dispersion, while J a b was varied until 
good agreement between simulation and experiment was 
achieved. 
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Following this procedure we determined that J a b = 
—6 ± 2meV. Figure [3{b) shows a slice through the sim- 
ulated data for J ab = — 6meV and J c = 12.2 mcV to 
give direct comparison with Fig. O^a). The distribu- 
tion of scattering within the plane is well reproduced by 
the model. Figure 0] shows constant-energy cuts through 
both real and simulated data sets along the line marked 
X in Fig. |21 at three different energies. The model does 
not include the variation of the background with energy, 
so a fiat background was fitted for each energy indepen- 
dently. In addition, the overall scattering amplitude had 
to be systematically reduced with increasing energy to fit 
the data satisfactorily. This reduction, which was nearly 
a factor of 2 over the energy range 6.5 meV to 14meV, is 
not predicted by the spin wave model. 

The measurements and calculations reported here re- 
veal that the magnetic correlations in Nao.7sCo02 are of 
a three-dimensional (3D) nature, despite its highly 2D 
physical properties. In fact, the inter-plane exchange 
constant J c is found to be roughly double the intra-plane 
constant J a b- The spin wave modes propagating along 
the c-axis are found to be sharp, indicating a well cor- 
related ground state. The in-plane modes exhibit some 
broadening, as indicated on Fig. 0] and reported previ- 
ously [Tol | . 

Comparisons have been made between Na x Co02 and 
other layered superconducting families, such as the cop- 
per oxides. The strong 2D nature of the cuprates is 
thought to be important for their superconductivity, and 
contrasts with the 3D magnetic interactions found here 
for Nao.75Co02. It is likely that the c-axis magnetic cou- 
pling is weakened in hydrated Na^Co02, due to the large 
inter-layer spacing, and it is tempting to speculate that 
this coupling actually inhibits superconductivity. This 
possibility is especially pertinent given recent evidence 
for the presence of H 3 + ions in the hydrated compound, 
which would make the doping level for superconductivity 
similar to that in Nao.75Co02 [l^ . 

The spin excitation spectrum observed here is not eas- 
ily reconciled with the usual picture of localized Co 4+ 
and Co 3+ ions carrying spins 5 = 1/2 and 5 = 0, re- 
spectively. If localized Co 4+ spins were distributed at 
random then a very broad magnetic excitation spectrum 
would be expected, unlike the sharp modes observed ex- 
perimentally. One possibility is that there is a phase 
separation into ferromagnetic in-plane clusters of Co 4+ 
ions in a matrix of non-magnetic Co 3+ . However, the 
Coulomb penalty would be considerable, and to obtain 
consistency with the observed sharp spin modes along the 
c-axis these clusters would have to be aligned vertically 
above each other over many layers. 

An alternative suggestion based on optical conductiv- 
ity data 01 i s that Nao.7sCo02 might have a stable 
Wigner crystal ground state, with Co 4+ spins on a tri- 
angular lattice of side 2a in a background of Co 3+ . This 
would double the period of the magnetic correlations 



within the layers and hence create magnetic zone centres 
at the M-points, e.g. (|00), of the Brillouin zone. Since 
no magnetic excitations are observed emerging from the 
M-points [10( we can rule out this type of charge order 
in Nao.75Co02, at least in the simplest case where the 
Co 3+ ions are non-magnetic. The Wigner crystal model 
can only be reconciled with the observed spin excitation 
spectrum if the Co 3+ ions carry a moment of similar size 
to the Co 4+ ions and all moments interact ferromagnet- 
ically. 

With these points in mind it might be that a more 
metallic picture should be sought, in which the Co atoms 
have no, or only a small disproportionation of charge. 
This would naturally produce a similar magnetic moment 
on each site, consistent with the model of an A-type anti- 
ferromagnet assumed here. In particular, a weakly itiner- 
ant ground state with strong spin fluctuations could rec- 
oncile the small ordered moment and low ordering tem- 
perature with the relatively large energy scale of the spin 
excitation spectrum, and account for the broadening and 
intensity loss of the spin modes with increasing energy. 
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